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Table 1: Ensembles (sets) of MILC configurations used with gauge coupling /3, size L? xT and sea masses 
(x tadpole parameter, mo) m 0l 9 and m 0s 9 . Column 3 is the lattice spacing values in units of r\ after 'smooth- 
ing' Q. Column 8 gives the number of configurations and time sources per configuration that we used for 
calculating correlators. On set 5 only half the number were used for light quarks. 

1. Introduction 

The B meson sector is a compelling target for lattice calculations for a variety of reasons. There 
are a variety of so-called "gold-plated" states — those states which are narrow and hadronically 
stable, as well as being experimentally accessible. In calculating these properties of states on the 
lattice there are no free parameters; m K , tuk, and mj calibrate the masses of the light, strange, 
charm and bottom quarks respectively, and T splittings and other meson masses calibrate the lattice 
spacing^, 

Precision B meson spectroscopy is a key ingredient in precision calculation of decay constants 
and form factors, ingredients in CKM matrix element determination and testing of the standard 
model. 

2. Simulation Methods 

We use five different ensembles of gauge configurations with 2+1 flavors of dynamical ASQ- 
TAD sea quarks, generated by the MILC collaboration. The ensembles, listed in Table [j], represent 
three lattice spacings, labeled very-coarse, coarse, and fine. 

On each configuration we generate and store random- wall HIS Q propagators for several source 
time slices for light, strange and charm quarks: 

S HISQ (x, ? o)=M->( ?0 ),', (2.1) 

where ri(to) x ' is a three-component complex unit vector of random numbers at each site of the 
source timeslice, to and zero elsewhere. 

The HISQ action uses an additional application of the fattening step of the ASQTAD formu- 
lation, reducing discretization errors to the extent that it is possible to simulate relativistic charm 
quarks on configurations of modest lattice spacing. 

Bottom quarks are too massive to simulate relativistically on these lattices. However within 
bound states, the b quark is generally slow enough (vb/c ~ 0.01 in B c ) to treat non-relativistically. 
The use of the NRQCD action for b quarks is a well-developed procedure. [Q, ||, ^] 
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We evolve the NRQCD propagator recursively: 



G ^ t + ^=\!--){}-fn) U ^{ l -fn) [ l --) G ^ ™ 
with 

fA( 2 )) 2 is - ~ is 

5H = - C ^ + C2 ^ {A - E - E - A) - C3 ^°- {AXE - EXA) 
g a 2 AW a (AP)) 2 

" C4 2M«° • B + C5 2Mfi ~ C6 T6^Y 2 • ( 23) 

The tilde expressions E and B are improved versions of the naive lattice chromo-electric and 
chromo-magnetic fields, E and B. We use the tree-level values of c,- = 1 for the constants. 

To double statistics, we evolve the NRQCD propagator both forward and backward across the 
lattice from the source timeslice. 

As we have used a random-wall source for the HISQ propagators, it is critical that we initialize 
the NRQCD b propagators with the same random-wall function T7(? )y as we use d for the HISQ 
propagators. This is slightly non-trivial in that the HISQ staggered fermions, and the random 
wall vector T](/o)y, have one Dirac component per site, while the NRQCD b quarks have two 
upper and/or two lower Dirac components. The trick is to undo the staggering transformation by 
multiplying the noise source T](?o)y at each site with the four-component staggering operator: 

a{x)=y*ttt 3 - (2-4) 

Furthermore, to isolate the meson ground-state, we smear the b propagator source with a 
Gaussian smearing function of varying radii r,-. Therefore, on timeslice to we initialize the NRQCD 
propagator as: 

Gf RQCD (x,? ) =£S(|*-/|;r ; )T],^o)n(y)r, (2.5) 

x' 

where T is an element of the Dirac algebra chosen to project out a desired meson state. 

At the sink end we must also multiply Q. (x) back into the HISQ propagator so that we can get 
a multi-Dirac-component object to trace with the NRQCD b propagator: 

G ms *(x,t) ah =g ms *(x,t )Mx,t) ab (2.6) 

Then our B meson correlator matrix is: 

C T {t-t,) ij = ^G m ^\xj)TS{\x-^\-r j )G^(x',t). (2.7) 

X 

3. Analysis 



We extract B meson energies from the matrix of correlators ( |2.7[ ) using a Bayesian factorizing 
fit to the form 

Cr(t-t ) i j = Y t °Vfijt*~ E * <f ~* ) + I b^b) M (-l)^e- E '^\ (3.1) 

k=\ k'=\ 
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where the second term fits the oscillating component inherent in staggered meson correlators. 

We look for high-confidence fits stable with respect to varying Af eX p, and t m [ n of the fit range. 
Where possible we simultaneously fit all the correlators coming from the same ensemble, to better 
account for correlated errors. 

In practice we fit B (light) and B s together in all cases except the fine ensemble (set 5). We fit a 
3x3 matrix of smeared correlators in all cases except for the B c fits on the very-coarse ensembles 
(sets 1 and 2). We always fit the pseudoscalar and vector states simultaneously. 

We are interested in the ground-state energies Eq and the ground-state of the oscillating parity- 
partner channel E' . A factor of 7075 relates the spin structure T of the direct channel with that of 
the parity partner channel, V. In this way a measured pseudoscalar correlator also contains a scalar 
meson correlator, and a vector correlator also contains an axial vector correlator at no extra cost. 

Because the relativistic relation between energy and mass does not hold for NRQCD b quarks, 
there is an unknown energy shift between the physical masses we are interested in and the the fitted 
energies. Instead we measure the splitting between the state of interest and a similar state with the 
same NRQCD quark content. 

We convert this splitting to physical units using r\ = 0.3133(23)fm Jj]], giving a 0.7% uncer- 
tainty in any measured splitting in our lattice calculation. Hence we can minimize the scale-setting 
error by choosing comparison states as close as possible to the state of interest. We are perfectly 
free to construct a fictitious comparison state which is a composite of real states, provided all 
components have well-known experimental and lattice measurement for calibration. 

We consider three methods to determine the B s and B c masses: 



Here E b ^, for example, refers to the spin-averaged lattice energy of bb states. In each equation 
we must apply the lattice scale a -1 (and its uncertainty) to the expression in the () latt only. 

Where the subtraction compares states with different electromagnetic charge structures we 
must estimate the adjustment necessary to account for electromagnetic effects. 

4. Results and discussion 

4.1 Pseudoscalar states 

In practice Method Q is the only one applicable to B s spectroscopy. We extract the lattice 
energies of the B s states from each of the ensembles, convert to physical masses via expression |. 

After the recent, more precise determination of r\ [|J], it has become apparent that both the s 
quark mass and the b quark mass were tuned too high. Method I for M Bs is particularly sensitive 
to the mistuned quarks. We have estimated the effect of the mistuned quark masses by substituting 
into Method | mesons with different valence masses. We estimate that for the very coarse, coarse 
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Figure 1: Lattice calculations for Mb s (left) Mg^ (right) and with energy shift subtracted and resulting 
masses extrapolated in a 2 to the physical point. For Mb, we correct the finite-a 2 points for b and s mass 
mistuning before extrapolation. For Mg ( , the continuum points are corrected upwards for electromagnetic 
effects by 4.5MeV and lMeV for Methods I and II respectively. Error bars on extrapolated points reflect 
total errors. 

and fine ensembles, the too-large strange mass pushes up Mg s by 7.5, 10 and 9MeV, respectively. 
The b mistunings bias Mg s up by 10.5, 13, and 15MeV on the same ensembles. We correct for 
these biases in the finite a calculations, and then extrapolate, estimating an additional systematic 
uncertainty of lOMeV on the extrapolated value, giving: 

M Bs = 5.341(4)(10)GeV, (4.1) 

with the first error being statistical and the second, dominant, error being the quark-tuning system- 
atic error. Figure |] (left) illustrates the extrapolation. 

For B c pseudoscalars we can use Methods || and [II As || is superior to Q and they are not 



linearly independent we do not also consider | here. We again extrapolate in a 2 to the continuum 
for each. 

We correct for the electromagnetic structure mismatch. Method || compares neutral bb and cc 
states with the charged B c state. We calculate that this mismatch causes an underestimate of M Bc by 
~ 4.5 ± 2MeV. In Method |l|, comparing similarly charged B c and D s introduces an underestimate 
of ~ 1 ± lMeV. 

After correcting the electromagnetic contribution we get: 

M Bc (II) = 6.279(2) (1) (5) (2)GeV (4.2) 
M Bc {UI) = 6.268(4)(6)(l)(l)GeV, (4.3) 

where the errors are (statistical)(ri )(NRQCD)(EM). The agreement between the two independent 
subtraction methods is a strong test of our control of systematics. Because the ()i att term is very 
small, both are quite insensitive to b and s tuning, and no further subtraction is necessary. Results 
from both methods are in excellent agreement with the PDG average of 6.277(6) GeV[||]. See 
Figure [I], right. 

The HISQ c quark seems to be the source of the strong discretization effects in Method II, 
which go as a s {v/c) r(am c ) . The c quark is more relativistic inside the B c than in a cc, so these 
errors do not cancel exactly, but should vanish in the continuum. 
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Figure 2: The extrapolation of the ratio R to the continuum for both B* and B*. 

4.2 Vector states 

For the B* states there is an obvious method of correcting for the NRQCD energy shift — 
compare to the nearby B pseudoscalar states to get the hyperfine splitting. The remaining compli- 
cation is that since the a • B term in the NRQCD action generates the B* — B splittings, radiative 
corrections to this term could generate a multiplicative correction to the splitting. (Recall we have 
used the tree-level C4.) We therefore use the hyperfine splitting of the B s system as calibration for 
that of the B c system and calculate: 

Er* — Er 

Rc= J : c F , (4.4) 

Eb* - E Bs 

which will cancel all of the NRQCD energy shifts, multiplicative corrections, and scale-setting er- 
ror. We extrapolate to a 2 = 0, multiply by the PDG average value of M B * —M Bs = 46.1(1.5)MeV [01, 
and add the experimental B c mass, giving us a prediction of the B* mass of Mb* = 6.330(7) (2) (6) 
GeV. As a check we also calculate /?/ with light quark B and B* states. A complete discussion of 
the B c hyperfine splitting calculation can be found in 

4.3 Scalar states 

As mentioned in Section 3, the oscillating component of the pseudoscalar correlators gives us 
the scalar states. 

We extract the E + — E - splittings directly in the simultaneous fits. Converting to physical 
units we again extrapolate in a 2 to the continuum and find: 

AM Bj (0 + -<r) =0.41(2)GeV (4.5) 

AM Bc (0 + -0~) = 0.44(7)GeV, (4.6) 

quoting statistical errors only. As this splitting is generated by the kinetic term, it should acquire 
no multiplicative renormalization, and most systematics should cancel. The + B c state lies about 
400MeV below the B+D threshold so it should be a narrow state. It is less clear whether the + B s 
state is below the B + K state, but in any case it should be close enough to the B + K that it should 
also be a narrow state. See Figure |3[ right and left, respectively. 
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Figure 3: The splitting between the scalar and pseudoscalar states for B s (left) and B c (right). Values are 
extrapolated to the continuum. Shown on each are the relevant hadronic thresholds. 

5. Conclusions 

We have shown preliminary results of precise lattice calculations of pseudoscalar masses in the 
B s and B c system, and of vector-pseudoscalar and scalar-pseudoscalar splittings. Our calculation 
of Mb s and Mb c agree within errors with experimental measurements of these states. Some work 
remains to fully understand the systematic errors and biases related to mistuning of quark masses. 

Our calculations of the B s and B c scalars and the B* vector constitute predictions of the masses 
of these states before experimental measurement. 

The precision and accuracy of these results reaffirms that the combination of HIS Q light quarks 
and NRQCD b quarks is a powerful lattice technique. Further work will complete the exploration 
of the lowest B,B S and B c states, and then apply these techniques to the calculation of form-factors 
and decay constants relevant to weak-matrix elements. 

References 

[1] C. T. H. Davies, E. Follana, I. D. Kendall, G. P. Lepage and C. McNeile, [arXiv:0910.1229 [hep-lat]]. 

[2] E. Follana et al. Phys. Rev. D 75, 054502 (2007) [arXiv:hep-lat/06 10092]. 

[3] E. Follana, C. T. H. Davies, G. P. Lepage and J. Shigemitsu Phys. Rev. Lett. 100, 062002 (2008) 
[arXiv:0706. 1726 [hep-lat]]. 

[4] A. Bazavov et al, arXiv:0903.3598. 

[5] B. A. Thacker and G. P. Lepage, Phys. Rev. D 43, 196 (1991). 

[6] G. P. Lepage et al, Phys. Rev. D 46, 4052 (1992) [arXiv:hep-lat/9205007]. 

[7] C. T. H. Davies et al, Phys. Rev. D50:6963 (1994) [arXiv:hep-lat/9406017]. 

[8] Particle Data Group, http://pdg.lbl.gov. 

[9] E. B. Gregory et al, arXiv:0909.4462 [hep-lat]. 



7 



